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Abstract 

We consider explicit models of dynamical supersymmetry breaking where dark matter is 
a 10 - 100 TeV strongly-interacting composite state carrying no standard model quantum 
numbers. These constructions are simple variants of well-known supersymmetry breaking 
mechanisms, augmented to allow for a large "flavor" symmetry. Dark matter is the lightest 
composite modulus charged under this symmetry and is a viable cold dark matter candidate 
with a thermal relic abundance. This is an attractive possibility in low-scale gauge-mediated 
scenarios where the gravitino is the lightest superparticle. A light i?-axion associated with 
supersymmetry breaking is present in these hidden sectors and serves as the portal between 
dark matter and the standard model. Such scenarios are relevant for present and future 
indirect detection experiments. 



1 Introduction 



The nature of the dark matter (DM) in the universe is a central open question in particle physics, 
astrophysics, and cosmology. A leading paradigm for the origin of DM is thermal freezeout of 
stable massive particle, which yields an appropriate relic abundance of cold dark matter [HIS]- 
The most commonly studied case is a weakly- interacting massive particle (WIMP), where DM 
has a mass of 100 GeV - 1 TeV and annihilates through weak couplings. However, there is a 
much larger range of DM masses and couplings consistent with thermal freezeout. 

In particular, strongly-coupled DM can have the right relic abundance if its mass is in the 
range 10 - 100 TeV. This can occur if DM is a stable state in a hidden sector which has small 
direct couplings to the standard model (SM) |31I11IS]. Indeed, there is independent motivation for 
this 10 - 100 TeV mass scale coming from supersymmetry (SUSY), specifically low-scale gauge 
mediation. An important ingredient for any model of low-scale mediation is a dynamical SUSY 
breaking (DSB) sector, and in the lowest-scale mediation models, the DSB scale is A ~ 10 - 100 
TeV. Here, we study the possibility that DM is a stable state in the DSB sector with no SM 
quantum numbers, with a mass determined by this dynamical scale. 

To set the stage, recall that gauge-mediated supersymmetry breaking (GMSB) [B] is an ap- 
pealing paradigm in its own right since it allows for natural suppression of flavor violation in 
the supersymmetric standard model (SSM) However, in low-scale GMSB, there is no good 
candidate for cold DM among the SSM fields. The lightest superpartner is almost always the 
gravitinojj which can only be a warm DM candidate if its mass m3/2 is about a keV§ So in the 
context of GMSB, a stable state in the DSB sector is compelling candidate for cold DM. 

Simply having strongly-coupled DM is insufficient for DM to be a thermal relic; there must be 
sufficiently strong interactions between the DSB and the SSM to maintain thermal equilibrium 
before freezeout. As advocated in Ref. |TU], strongly-coupled DM can couple to the SSM via 
a light mediator field, dubbed a "portal". Since DSB is almost always associated with sponta- 
neous breakdown of an R symmetry [11], a natural candidate for this portal field is the light 
i?-axion. Moreover, if DM is only quasi-stable with a lifetime r ~ 10^^ sec, then DM can decay 
via this i?-axion [lOj, potentially explaining the electron/positron excesses seen in cosmic ray 
observations [El [I3l [H [15] . 

In this paper, we study DM models based on different DSB mechanisms, and discuss the 
necessary conditions for the existence of a viable DM candidate. In addition to having a portal 
field like the i?-axion, DM in the DSB sector requires a global symmetry preserved by the strong 
dynamics, such that DM can be the lightest state charged under this symmetry. While the 
simplest DSB models frequently considered in the literature do not satisfy these conditions (e.g. 

^There are notable exceptions, see Refs. [8]. 

^In addition, structure formation typically requires the gravitino mass to be lighter than 0(10 eV) [S]. 
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the 3-2 or 4-1 models [6l [161 [U])) 'we find there are straightforward generahzations in which 
composite DSB states can be viable DM candidates. 

Concretely, we consider two explicit models with viable DM based on two popular DSB mech- 
anisms: 

• SUSY breaking from a dynamically generated superpotential P, [161 IE] ; 

• SUSY breaking from a quantum moduli space (QMS) |T8]. 

We also briefiy mention DM in the context of metastable SUSY breaking [19J. Since we will be 
considering these DSB sectors in the strongly-coupled limit, we will only be able to study the 
symmetry structures of these sectors and estimate their DM properties. However, finding DSB 
sectors that satisfy the necessary conditions for composite DM is already interesting, and we 
expect that these requirements will open new directions in DSB sector model building. 

We will show that under reasonable assumptions about the DSB sector, DM can be cosmo- 
logically stable and have the desired thermal relic abundance. While not the main focus of this 
work, we will also show that such models may be relevant for indirect detection experiments, in 
particular the possibility of DM decay being the source of several recently observed excesses in 
PAMELA [12j, FERMI [13] and H.E.S.S. [H [B] (see Refs. [101 1201 [211 1221 [23] ) . This occurs if 
the symmetries in the DSB sector that stabilize DM are violated by physics at the unification or 
Planck scale, and our explicit constructions naturally accommodate this possibility. 

This paper is organized as follows. In Sec. [21 we discuss the desired properties of DM from 
the DSB sector and briefiy mention alternative scenarios. We then construct two explicit models 
in Sec. [3] and Sec. [H showing that two popular DSB mechanisms allow the necessary symmetry 
structures to have viable DM. In Sec. [5l we briefiy discuss how such models can fit into the GMSB 
framework by introducing messenger fields. We conclude in Sec. [B 

2 Viable dark matter from DSB 
2.1 Minimal requirements 

Our starting point is a DSB sector with fields and interactions characterized by a scale A ~ (10 
- 100 TeV). This scale is determined by the strong dynamics that triggers SUSY breaking. In 
order for the DSB sector to have a DM candidate, it must possess at least one (accidental) global 
symmetry which is left unbroken by the strong dynamics. The lightest state charged under this 
global symmetry is a stable DM candidate, and the natural size for its mass tudu is the SUSY 
breaking scale A. We briefiy mention other possibilities in Sec. 12.41 

The relic abundance of such a composite is determined from its annihilation cross section into 
other lighter strongly-interacting particles not carrying the conserved quantum number, which 
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are in turn kept in thermal equilibrium with SM degrees of freedom through some kind of portaljfl 
Assuming 2—7-2 annihilation, the thermally averaged cross section is 



57r m 



DM 



where k denotes the relevant coupling of the hidden sector. As recently emphasized in Ref. [TU] , 
for typical couplings of a strongly-coupled sector k ~ (A/ivr — 47r) and mDM ~ 0(10 — 100 TeV), 
we obtain the desired relic abundanceo With SUSY breaking (and messenger fields) at this scale 
a/F ~ C(10 — 100 TeV), gauge mediation generates MSSM gaugino and scalar masses of order 
(^Vl67r2)/F ^ 0(100 GeV - ITeV). The gravitino mass is mg/a = F/VSMpi ^ (0.1 - 10 eV), 
where Mpi ~ 10^^ GeV is the reduced Planck scale, and the resulting gravitino abundance is small 
and consistent with structure formation |9]. 

A simple way to keep the DSB sector in thermal equilibrium with the SSM is to assume that 
both sectors are charged under the same (approximate) global symmetry which is subsequently 
broken spontaneously [lOj. In SUSY breaking models, R symmetry is a natural candidate, and 
we will focus on that case throughout. Generic DSB models will spontaneously break this R 
symmetry [11], and as long as there are no large sources of explicit R violation, the resulting 
i?-axion will be light. The i?-axion is then coupled to both the DSB and SSM sectors, and 
will keep the two sectors in thermal equilibrium. Intriguingly, supergravity yields an irreducible 
contribution to the i?-axion mass, since to cancel the cosmological constant, the superpotential 
must have a R symmetry breaking constant term FM-p\ . After taking into account of possible 
additional explicit R symmetry breaking effects, the i?-axion mass is then in a phenomenologically 
interesting range ma ~ 0{1 MeV — 10 GeV), with implications for collider experiments and dark 
matter indirect detection [TP] . 

To summarize, the minimal symmetry requirements to have strongly-coupled thermal relic 
DM in the DSB sector are: 

1. An (approximate) unbroken global symmetry under which DM is charged. 

2. A spontaneously broken (approximate) global symmetry under which both the SSM and 
DSB sectors are charged. The natural candidate is an R symmetry, resulting in a i?-axion. 

Note that this second requirement of an "axion portal" is not unique. For example, one 
could maintain thermal equilibrium between the DSM and SSM sectors using the vector U{1) 
portal [27], which we will briefly mention in Sec. 



■^In GMSB models, thermal equilibrium might also be maintained through couplings to the messenger sector. 
^The right relic abundance could also be maintained if the DM mass is lower than the weak scale, and if the 
relevant couplings is smaller, see Ref. [24] , 
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2.2 Strong coupling 

For the large values of k necessary to get the right relic abundance, one might worry about 
maintaining theoretical control over the DSB dynamics. In general, one is not able to calculate 
the Kahler potential in a strongly-coupled SUSY theory. Thus, the details of the low-energy 
theory below the strong dynamics scale — such as the masses of low-lying states and the proper 
identification of unbroken global symmetries — are strictly speaking unknown. Our strategy is to 
identify DSB models which satisfy the necessary conditions to allow a DM candidate, and assume 
that the symmetry structure at small values of k is preserved at strong coupling. 

Indeed, it is already non-trivial to identify explicit weakly-coupled DSB models that possess 
the necessary symmetries to have stable DM. While the above requirements in Sec. 12.11 were 
already outlined in Ref. [10], that work could not identify an explicit single-scale DSB model with 
viable DM. In particular, Ref. [10] was only able to show that a dynamical (but SUSY-preserving) 
sector could be consistently coupled to an i?-breaking and SUSY-breaking O'Raifeartaigh model. 
In the present work, we will show two examples of weakly-coupling DSB sectors with the right 
structure to contain DM, and then extrapolate to strong coupling. 

A key assumption we will make in our studies is that the relevant degrees of freedom for 
describing DM are gauge-invariant moduli fields. However, since we are considering the limit 
where moM — A, in principle other composite degrees of freedom might be the true DM modes. 
That said, since there is always a parametric limit where the moduli can be made lighter than 
the other composites, we think this assumption is justified. 

2.3 Dark matter detection 

Though not the main focus of this paper, we will show that DSB DM models may be relevant for 
DM indirect detection experiments. In particular, recent cosmic ray excesses [121 lISl ttH [15] might 
be explained by heavy DM decay to the i?-axion [10]. This can occur if the global symmetries 
protecting DM are only accidental at low energies, and are not respected by physics at high scales 

such as the unification or Planck scale [22] . In the low energy effective theory, these breakings 
are encoded in high dimension operators suppressed by M,,, which can induce DM decay, and we 
will see that DSB models naturally accommodate such operators. We emphasize that it is possible 
that such operators are absent in the DSB DM models considered here. In this case, the stable 
state in the model still provides viable DM candidate. Although it is not possible to explain the 
recent cosmic ray excesses, the annihilation of dark matter can still give indirect detection signals 
which may be probed by future observations. 

The DM lifetime which can potentially explain recent cosmic ray excesses is of order 10^^ 
sec [22]. For 1 - 100 TeV DM, dimension-six symmetry breaking operators suppressed by the 
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unification or Planck scale could explain the anomalous astrophysical signals [22]. For i?-axion 
decay constants in the range fa ~ 0{1 — 100 TeV), the most natural region avoiding all the 
constraints from astrophysics and rare meson decays are > 2m^ [lOj. For the range ^ 
(200 MeV — 10 GeV), the i?-axion decays mostly into either t^t~ , or tt+tt^tt^ and provide 

a good fit for the electron/position signals in cosmic rays [TO] . 

In the scenario considered in this paper, it is difficult for direct detection experiments to 
discover DM if it is not charged under SM gauge groups. The SM-singlet DM can only interact 
with nuclear matter via i?-axions or loop diagrams mediated by messenger fields, and the resulting 
cross sections are unobservable with present techniques. 

2.4 Alternative constructions 

In this paper, we will focus on the possibility that DM is a 10 - 100 TeV state in the DSB sector 
with no SM quantum numbers. Here, we briefly mention some alternative possibilities for DM in 
GMSB-like scenarios. 

First, the messenger sector — either a separate sector or part of the DSB sector — has previously 
been proposed as a possible source for a cold DM candidate in GMSB models [3l [5l [281 129] . 
The lightest messenger particle could be stable provided the existence of a "messenger number" 
invariance. Assuming that they only annihilate into particles of the SSM through SM gauge 
interactions, they should be lighter than about 5 TeV in order not to overdose the Universe. 
Direct detection experiments have already ruled out such a WIMP [30] unless the elastic scattering 
of DM off nucleons is forbidden [2B]. Another possibility is to have a strongly-coupled composite 
heavy messenger at the scale of 10 - 100 TeV as DM [21]. As we mention in Sec. 0, even if DM 
dominantly arises from the DSB sector, the messenger sector often has stable states of its own, 
which must be taken into account. 

In our constructions, we focus on the DSB sector, and we will be considering the case that 
the mass of DM mass is comparable to the dynamical scale. However, it is also possible for 
the DSB sector to contain stable states which are parametrically lighter than the DSB scale, 
e.g., pseudo-moduli or pseudo Nambu-Goldstone bosons (pNGB). They can potentially be good 
cold DM candidates, as has been recently explored in the literature [^ [52] 155] . From the basic 
requirements for DM such as the correct relic abundance, though, it seems necessary to extend 
beyond the basic ingredients of the DSB sector in order for such scenarios to be viable. 

For example, a weakly-coupled Intriligator-Thomas-Izawa-Yanagida (ITIY) model with scalar 
mesons as DM has been studied in [31] . The scalar mesons are pNGBs with suppressed couplings 
A ~ C(10~^ — 10""^) and masses around TeV. To achieve the right amount of relic abundance, a 
"resonance mechanism" was invoked by tuning the mass of the pseudo-modulus S to be very close 
to twice the meson DM mass. Once the proper relic abundance is achieved, possible DM decays 
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in Ref. [31] were induced by a holomorphic Kahler potential terms K ~ AMj or equivalently 
dimension five superpotential terms after taking into account of supergravity, similar to what we 
will do in Sec. 14.31 

Both pseudo-moduli and pNGBs could be DM candidates in models with metastable SUSY 
breaking [321 [33]. It has recently been shown by Intriligator, Seiberg and Shih (ISS) that deformed 
SQCD theories with weakly-coupled magnetic duals have a meta-stable SUSY breaking vacuum 
near the origin of field space [H]. The original ISS model has massless fermionic moduli fields and 
some pNGBs from flavor symmetry breaking (baryons or dual quarks). Both are charged under 
flavor symmetries (baryons are further charged under a U{1)b) and the lightest states are stable. 
The annihilation of the moduli to pNGBs or vice versa will be suppressed by couplings much 
smaller than and the resulting relic density tends to overdose the universe. Thus to have 

a viable DM in ISS models, one has to make the fermionic moduli massive and have additional 
interactions for DM to annihilate more efficiently. 

We will not try to build a detailed DM model based on the ISS model but just point out 
possible solutions. First, to make the fermionic modulus massive, one must break an accidental 
R symmetry in the ISS model. This can be understood in analogy to the gaugino masses in 
the gauge mediation, since both fermion mass terms break a continuous R symmetry. Thus 
both fermions can acquire masses at one loop after R symmetry is broken. R symmetry can be 
broken either explicitly by additional fields and operators or spontaneously through the "inverted 
hierarchy" mechanism [3^ . 

Second, to have large enough DM annihilation cross section in the ISS model, one could gauge 
some unbroken flavor symmetries. For instance, in the Nf = N^+l ISS model, we can adjust the 
parameters such that the baryon is the lightest field charged under the global symmetries. Then, 
we can gauge the U{1)b which can kinematically mix with the SM hypercharge U{1)y- The DM 
sector then communicates to the SM through the U{1) vector portal [23132], and it is possible to 
have the right DM relic density. Another model is constructed in Ref. [35] where the unbroken 
flavor symmetry is gauged and identified with the SM gauge group. The fermionic pseudo-moduli 
then annihilates into the SM fields through SM gauge interactions. 

3 Dark matter from a dynamical superpotential 

In this section, we show that the requirements from Sec. 12.11 to achieve DM in the DSB sector 
can be met in models with a non-perturbatively generated superpotential. Many early models 
of DSB were based on such constructions, e.g., the 3-2 model [I6] and the 4-1 model [6l [T7] . 
In these models, a non-perturbative superpotential lifts the origin of field space while an added 
classical superpotential constrains some fields to stay near the origin. SUSY breaking can be 
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established by the Affleck-Dine-Seiberg (ADS) criterion fl6\ [36] : the classical scalar potential 
lifts all non-compact flat directions and there exists a spontaneously broken global symmetry (an 
R symmetry in this case). 

The existence of a spontaneously broken R symmetry meets one of the criteria for DM in the 
DSB sector. However, the 3-2 and 4-1 models do not have additional unbroken global symmetries 
that would lead to stable DM. As we will see, by extending the symmetry structure of these models 
one can get viable composite DM. Our construction requires a large unbroken global symmetry, 
so in this sense, it is similar in spirit to "single sector" models where SSM degrees of freedom are 
composites from a DSB sector [37] . 

3.1 General considerations 

Recall that the 3-2 and 4-1 models are the simplest examples of two infinite classes of DSB models 
based on SU{n) x SU{2) x f/(l) (n odd) and SU{n) x U{1) {n even) gauge groups. They can 
be constructed using the discarded generator method [51 EH]- For example, SU{n) x U{1) {n 
even) models can be deduced from the non-calculable theory with gauge group SU{n + 1), an 
antisymmetric tensor, and n — 3 antifundamentals. After discarding those generators of SU{n + l) 
which do not lie in the subgroup SU{n) x f/(l), the original chiral fields decompose as 



where A is an antisymmetric tensor, F{F) is the (anti)fundamental, S* is a singlet of SU{n), and 
the subscripts denote the U{1) charges. 

When n > 6, these models have an SU {n — 3) global symmetry under which F and S are 
charged. At the classical level, one can add to this model a superpotential which lifts all flat 
directions 



For suitable choices of the A and rj matrices, the global SU {n — 3) (or some subgroup) will be 
preserved. If the symmetry is not broken spontaneously after SUSY is broken, this would signal 
the existence of stable degrees of freedom to be DM candidates. The same argument can be 
applied to the other models constructed in this way. 

The minimal models, i.e. 3-2 and 4-1, possess no stable degrees of freedom as there is no 
global symmetry left after SUSY breaking. Thus, in order to incorporate a DM candidate, the 
field content needs to be enlarged to allow for a larger global symmetry which remains unbroken 
after SUSY breaking. 



A2 + Fi-n + {n- 3)F_i + (n - 3)5, 



(2) 



W = XijAF'F^ + r]ijFF'S^. 



(3) 
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3.2 An explicit SU{Q) x U{1) model 

The simplest example in SU{n) x U{1) {n even) with a non-Abelian global symmetry is the 
SU{6) X U{1) model, and this model does have the right symmetry structure to have stable DM. 
The microscopic field content of the SU{Q) x f/(l) model is: 





SU{6) 


m 


SU{3) 


Uil)R 




15 


2 


1 


-4 


pa 


6 


-5 


1 


3 


K 


6 


-1 


3 


3 


Si 


1 


6 


3 


-4 



where the Greek scripts are the gauge indices and the Roman scripts i {— 1,2,3) are the global 
SU (3) indices. At the classical level, one can add a superpotential: 

Wa = Xer^sA-f'F^F'^ + ViF'-iF^S, + F^S^) + VsF^F^Ss, (4) 

which lifts all the D-fiat directions as we will show later. This classical superpotential explicitly 
breaks the global SU{3) to SU{2). Notice that although we choose to preserve SU{2), a global 
U{iy or U{1) would be sufficient to satisfy the minimal requirement for DM inside the DSB 
sector. 

For a range of parameters of the model A, r/ <^ 1, the fields obtain vacuum expectation values 
(VEVs) much above the dynamical scale A of the SU{6) theory. The theory is then weakly- 
coupled and one could analyze SUSY breaking in terms of the microscopic variables since the 
Kahler potential of the light degrees of freedom is nearly canonical. Instead, we focus on the 
strongly-coupled regime where the scale of SUSY breaking is around the dynamical scale A. In 
the effective theory below the scale A, the relevant degrees of freedom are the gauge invariants, 
i.e. moduli, built out of the microscopic fields. Note that U{1) is infrared free and its gauge 
coupling is weak at the scale A. 

The SU (6) invariant moduli fields are 

S„ M' = F^Fl H, = e,,kA''^FiFl 5 = Pf A = e„,,,,3,,«3«,A"i"^A"^°^A"^-«. (5) 

They have charges 15(6), M(— 6), H{0), -8(6) under the gauged U{1). In general, one expects non- 
perturbative SU{6) dynamics to cause M or _B or both to gets VEVs, breaking the U{1) gauge 
symmetry. The relevant degrees of freedom in the low energy theory are the U (1) invariants built 
out of the above SU (6) moduli fields: 

Xi^SjM\ Hi, Y'^M'B. (6) 

They are not all independent and satisfy the classical constraints: 

^prk _ -^kyi ^ 0. (7) 
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It can be shown explicitly by global symmetries that the constraint is not modified quantum 
mechanically|§ Notice that the 9 constraints in Eq. ([7]) are not all independent. For each fixed 
j, only two of them are independent. In total, there are 6 independent constraints which allow 
us to eliminate 6 degrees of freedom from the 15 moduli fields listed in Eq. (Q. We choose the 
independent 9 degrees of freedom to be 

X|:(2 + l)2, Hr.{2 + 1)2, r^(2 + l)_6, (8) 

where we have listed the quantum numbers of the moduli under the unbroken global SU{2) x 
U{1)r. Note that we have 10 independent SU{6) moduh in Eq. ([5]) before gauging f/(l). After the 
gauged f/(l) eats a multiplet, there should be only 9 independent fields left, which is consistent 
with our explicit construction and counting above. 

In terms of the gauge invariants, we express the classical superpotential as 

W,i = XH, + r],{Xl + X^)+r]sXl 



= Aif3 + r^i(^^^^^^^^)+r^3X|, (9) 

where in the second line, we have expressed everything in terms of the independent moduli fields. 
Now consider the microscopic equation of motion (i.e. in terms of A, F, F, and S) corresponding 
to this classical superpotential. The equations dW/dS sets X, Y to zero if we multiply by S, B. 
Similarly dW/dF sets H to zero if we multiply by F. Thus at the classical level, our theory has 
no fiat directions, satisfying one of the ADS criteria. 
The dynamically generated superpotential is [17] 

W.yn = (10) 



The exact superpotential is then a combination of Eqs. Q and ffTOj) . The equation of motion for 
X| cannot be satisfied, and thus SUSY is broken. In the regime near the origin of the moduli 
space, the Kahler potential is smooth in terms of the composite fields, e.g., K^s ~ + l-f^P- 
Minimizing the potential shows that there is a runaway direction with ^ oo. But if (Y) is 
large compared to A, we should treat the quarks as elementary degrees of freedom and the Kahler 
potential is smooth in A, F, F. It can be shown that the potential rises when {A, F,F) ^ A [B]. 
We see that there is no supersymmetric vacuum for either large or small values on the moduli 
space, and so there could be at least a local SUSY breaking minimum for {Y, H) ~ A where the 
U{1)r symmetry is also broken. 

Because the theory is strongly coupled, the Kahler potential for the moduli is in general 
unknown. So while we have established that SUSY and U{1)r must be broken, we cannot say for 



'That is, XY has i?-charge —4 so it is inconsistent to have a zero i?-charge instanton term in the constraint. 
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certain what the fields VEVs are. We will assume that there is a (local) minimum at the branch 
where the global SU (2) remains unbroken, which means: 

(ffg, X|, Y^) ^ 0, other VEVs = 0. (11) 

This is a reasonable assumption, as points of maximal symmetry are generically the stationary 
points of the energy. But one should bear in mind that this could be modified due to our lack of 
knowledge about strongly-coupled theories. 

3.3 Spectrum and annihilation 

The stable DM candidate is the lightest component(s) of Xf,Y^,Hi with i = 1,2. Since we 
are in a regime where the theory is strongly coupled, it is not possible to directly calculate the 
annihilation cross section for the DM candidate. However, we will show that the annihilation 
cross section is generically of the right size to get the desired relic abundance. 

First note that the fields that get VEVs in Eq. f lTTj) are indeed unstable. The R symmetry 
is broken spontaneously, thus the model has a light i?-axion with the decay constant fa of order 
A. The couplings of the modulus to the i?-axion are fixed by sjTiimetries and can be read 
off from the modulus kinetic terms [1^. For the scalar components, the kinetic term is £ = 
+ 2r$(($) + $')(9^a)//ap, where $' = $ — ($) in which $ represents a the scalar field 
with non-zero VEV, such as H^, X|, and Y^. r$ is an order one coefficient depending on the 
i?-charge of that field. Assuming that all the VEVs are real, this yields a coupling of Re($) to 
two i?-axions, and coupling between Re($), Im($) and an i?-axion. Therefore, Re($) decays 
promptly into two i?-axions, while Im($) decays to three R-axions via an off-shell Re($). SUSY 
breaking leads to an coupling of the fermionic component of $ to its scalar component and a 
gravitino. This allows the fermions to decay into a gravitino and an i?-axion. 

The massive U{1) vector multiplet is also unstable. After being higgsed, the f/(l) gauge 
boson and gaugino will obtain a mass of order ~ giA/{47i) where gi is the U{1) coupling. 

As Qi ^ Att, they are lighter than most of the composites with mass around A except for the 
i?-axion. After the U{1) is broken, gauge-variant operators like K = giVl^l"^ can be generated 
with V being the vector multiplet. Expressing the operator in components, one can see that it 
causes the gauge boson to decay to one i?-axion and the off-shell unstable scalar $, which will 
decay further to i?-axions. The gaugino will decay to the off-shell fermionic and scalar $, both 
of which will cascade to i?-axions. 

Contributions to the DM spectrum and annihilation cross section come from two sources: the 
superpotential and the Kahler potential. A detailed spectrum is impossible to obtain but we can 
still estimate the sizes of the masses and interaction strength. After canonically normalizing the 
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composite fields, the superpotential becomes 

W = ~XA^H, + f^,A^ (^^^^^^^) + mA^Xl + (12) 

where the coefficients A, f/, a absorb unknown factors from normahzation. Using naive dimensional 
analysis (NDA) |39], the sizes of the coefficients in the superpotential and the field VEVs are 

{Fh,x,y) ~ 9^', {Hs,Xl,Y') ^gA, (13) 

where g ~ (47r)~^. We will take {X,1]) to be of order one and thus {X,fi) of order g ~ (47r)^^. 
Expanding around the minimum of the potential, the superpotential gives mass terms for the 
composites as 

W ^ A{Y'H, + Y^Hs + XfY'), i = 1,2, (14) 

where we neglected order one coefficients. These leads to supersymmetric masses of order A ~ 
(10 — 100 TeV) for nearly all the fields except for X| and a linear combination of Xf and Hi. But 
these fields will receive higher order Kahler corrections we will discuss below. The interaction 
terms from the superpotential are 

W ^ g-\Y'Hi{Y'' + H^) + X^Y'Y^), i = l,2, (15) 

where again we neglected order one coefficients and the interaction strength is of order ~ 47r. 

As the theory is strongly-coupled, control over the Kahler potential is lost. In general, the 
composites will feel 0{1) SUSY breaking effects depending on the details of the Kahler poten- 
tial. For instance, there could be higher-order Kahler terms compatible with all the unbroken 
symmetries as 

\Y\\ \X\\ irpli/p, (16) 

which are suppressed by the QCD scale A. These operators will give non-supersymmetric masses 
to all composites of order A and interaction between them with couplings of order 47r. For instance, 
will give masses to the scalar component of X| and {XfX^'^Y will lead to Majorana mass 
terms for Xf. 

DM could annihilate into the light vector multiplets. The annihilation cross section is sup- 
pressed by the weak couplings g\ except for the annihilation into longitudinal modes of the gauge 
multiplet V . In that case, the gauge coupling suppression is compensated by the polarization 
vector product (A^/m^)^ ~ {AiiY / g\ where we take my ~ (7iA/(47r). In addition, if the unstable 
fields (e.g. Y^) are lighter than the DM particles (e.g. Y^), the annihilation will proceed through 
diagrams like — ). Y^Y^. With strong couplings from NDA estimates, one expects to obtain 
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the correct relic abundance from such processes. Notice that we have not included multiplicity 
factors such as the color factor N in our NDA analysis, which will in general decrease the cou- 
plingsj^ Thus we do not prefer to work at very large A^. On the other hand, the thermal freezeout 
calculation will be affected by various co-annihilation channels if there exist states of comparable 
masses [IQ], which may partially compensate for large suppressions. Also, even if the unstable 
particles are heavier than the stable DM particles, thermal freezeout can still proceed, as long as 
the mass splittings are not too large [lU] . 

3.4 Dark matter decays 

An intriguing feature of these heavy composite DM models is that the global flavor symme- 
try is generically only an accidental symmetry at low energies. In particular, higher dimension 
operators — presumably encoding physics at the unification or Plank scale — can violate the flavor 
symmetries (and the f/(l)R symmetry) while preserving the gauge interactions. This allows for 
the possibility of DM decays, but on cosmological timescales. In our scenario, the decay occurs 
through the light i?-axion as in Ref. [10], so the final states from i?-axion decays can be mainly 
leptons, explaining the electron/positron excess observed in PAMELA, H.E.S.S., and FERMI. 
The required lifetime of order 10^^ sec is obtained if the decay is caused by a dimension six op- 
erator suppressed by the unification or Planck scale. If these dimension six operators break R 
symmetry, they will also contribute to the i?-axion mass. However, as the contribution is tiny 
compared to the irreducible supergravity effect, the i?-axion remains light. 

The global SU (3) symmetry was explicitly broken to SU{2) by the classical superpotential in 
Eq. (jlj). Assuming this SU{2) is preserved by the strong dynamics, the leading S*?/ (2)- violating 
gauge-invariant operators in this model are: 

W = ^FF^FfA, K=^ie,,sAF^F-' + S,FF'), (17) 

where is the high scale. We show the flavor indices explicitly and take all coefficients to be of 
order one. Below the scale A, these operators are matched into 

W=^c,Y\ K = ^^c,{H, + Xf), (18) 

where i = 1,2. The holomorphic Kahler terms can be translated into superpotential terms as 
W = Ci{A'^/M^Mpi){W + Xf). After including these two terms, the VEVs of the DM fields will 
be shifted away from zero to: 

{H,,Y\Xf) = c,^, (19) 



^Since it is the loop expansion factor k'^N that saturates the unitary bound (47r)^, the couphng k in Eq. ([T]) is 
suppressed by 
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where the Cj are coefficients of order 0(l/47r) by NDA. 

The tiny VEVs of the composites iJ, Y, X will induce them to decay into the i?-axion. If the 
DM is the scalar component of the composite, we have H,Y,X^ aa. For fermionic DM, one a 
will be replaced by a gravitino. For the scalar DM decay, the lifetime was estimated as [10] 

fn 27 /V47r\V \VlOTeV\^ , , 

-oM = 8.-^«10-=ec(^^j (lolMSv) (— ) (20) 

where w is a typical field VEV as in Eq. ( TT9l) with c a general coefficient of size q. We have also 
used ~ 47?/^ ^ A for the final parametrization. 

Thus, we have achieved an explicit DSB model where DM is stabilized by an approximate 
symmetry of the DSB sector. DM annihilates through unstable states which decay to i?-axions. In 
addition, DM itself can decay on cosmological scales through small symmetry breaking operators. 
In Sec. ini we explain how this DSB sector can be coupled to a messenger sector to have a complete 
model of gauge mediation. 

4 Dark matter from a quantum moduli space 

In addition to the chiral models discussed in Sec. |3l there is also a class of vector-like DSB models 
based on theories with a QMS, e.g. the ITIY model [18]. In these models, SUSY QCD with a QMS 
is coupled to some gauge singlet fields S with a tree-level superpotential SM, where M represents 
some composite fields from the SQCD sector. The equations of motion with respect to the singlets 
force the SQCD degrees of freedom to be at origin, which is not on the QMS. Thus SUSY is broken. 
These models usually accommodate large fiavor symmetries and some stable fields, which satisfy 
one of the minimum requirements of DM inside DSB models. The main complication in ITIY-like 
models is to establish the spontaneous breaking of an R symmetry, which is essential for both 
the generation of the MSSM gaugino masses and for DM communicating to the SSM through the 
i?-axion portal. 

Below we will study the simplest model of this class, the ITIY model. Such models have 
previously been studied in the context of DM in Ref. [31], albeit in the weakly coupled regime. 
Here we will focus on the strongly coupled regime. 

4.1 The ITIY model 

The ITIY model contains an SU{2) gauge symmetry with four fundamentals Q, as well as six 
singlet fields Sa- The fundamentals enjoy an S'f/(4) = S0{6) global symmetry. The field content, 
both in terms of microscopic fields and SU{2) invariant "mesons" is 
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where a = 0,1,2, ...5, and we have indicated the non-anomalous R symmetry. To this model, 
we add a classical superpotential of the form Wd = X'^SiMj. To satisfy the quantum modified 
constraint, at least one of the Ma will have to acquire a VEV, which spontaneously breaks the 
global symmetry 5*0(6) — t- 5*0(5). In order to remove the resulting massless Goldstone mode, 
the superpotential must contain explicit global symmetry breaking. The minimal choice, which 
we will make for our model, is 

PFd = Ao5oMo + Ai5iMi, i = l,...h. (21) 

It explicitly breaks the global 50(6) symmetry down to 50(5), and we assume that the coefficients 
A are order 1. Additional global symmetries can be explicitly broken by choosing a more general 
To achieve stable DM, it is important to assume that the potential at least preserves an 
accidental global 50(2) flavor symmetry. 

The equations of motion for the 5 fields set all meson VEVs to be zero, which contradict the 
quantum modified constraint 

MqMo + ^i^i = (22) 
i 

Thus SUSY is broken. Analogous to Sec. 13. 2[ we assume that the vacuum preserves the maximal 
50(5) global symmetry, such that of the meson, only Mq obtains a VEV (Mq) = A. Therefore, 
the fields M' = Mi and 5' = 5^ for z = 1, ...5 are stable DM candidates since they are charged 
under the (unbroken) 50(5) global symmetry. 

Note that Mq carries no i?-charge, so unlike the 6-1 model in Sec. El there is no guarantee 
that the i?-symmetry is broken. Below, we will make the assumption that Sq does get a VEV, in 
order to have the desired phenomenology. 

4.2 Spectrum and annihilation 

One can solve the constraint in Eq. (l22l) by writing Mq = (A^ — M'^)^/^. The resulting effective 
superpotential is 

Weff = Ao5o(A2 - M'2)i/2 + Ai5'M', 

= AoA25o- A[,5oM'V2 + AiA5'M' + 0(M'^), (23) 

where again S' = Si, M' = Mi for i = 1, ...5. In the second line, we have canonically normalized 
the composite fields and the coefficients A absorb unknown factors from the (non-canonical) 
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Kahler potential. In this form, is clear that the ITIY model is a type of O'Raifertaigh model. By 
NDA, we have 

Ao ~ ^, 9~\ Ai ~ 1, (24) 

which is obtained after we take the A, to be order one and g ~ (47r)^^. 

After SUSY breaking, there is one linear combination of the singlets S whose VEV is un- 
determined assuming the minimal quadratic Kahler potential for 5*. This is the tree-level flat 
direction of any O'Raifertaigh model. Previous analysis [H] has shown there is a minimum at 
(S) = by calculating the one-loop contributions from the light particles. Such an analysis is 
valid in the weakly-coupled regime, but for strongly coupled DM this conclusion may be modified. 
For example, when (5*) ~ A, the one-loop analysis breaks down as all loops contribute at the 
same order. Thus one cannot exclude the possibility that there is a (local) R symmetry breaking 
vacuum in the strongly-coupled region {S) ~ A, as assumed in Ref. [HI |12]. We assume that 
such a vacuum exists with (Sq) = ~ A while (5") = 0, such that the model still preserves the 
original SO (5) global symmetry. As in Sec. [31 the unstable field 5*0 will then decay to i?-axions. 

Expanding around such a vacuum, the effective superpotential is 

Weff = AoA^^o - Ao(w + ^o)M'V2 + AiA^'M', (25) 

where 5*0 = So — v. We see that S' and M' get a vector-like mass of order A. The field 5*0 will 
receive masses from higher order Kahler corrections such as |S'o|^. Assuming that 5*0 is lighter than 
the lightest component of M', the annihilation channel could be M'M' — )■ 5*0 5*0. The interaction 
terms generated from the superpotential such as V" ~ Ao^l^oPlM'p + ■ ■ ■ lead to an annihilation 
cross section of the order A[f/A^ ~ (47r)'^/(100 TeV)^ which has the right parametrics to yield the 
correct dark matter relic abundance. 

Thus, as long as we assume that a strongly coupled ITIY-type model dynamically breaks an 
R symmetry, then it has the right phenomenology to produce DM in the DSB sector. 

4.3 Dark matter decays 

As in Sec. 13. 4[ DM could decay on cosmological timescales via flavor- violating operators sup- 
pressed by a high scale. The leading flavor-violating operators are the Kahler potential terms 

K = ^.vMQVlQ', (26) 

with arbitrary flavor structure for r/"^. Below the confining scale A, this operator is mapped onto 

K = ^c^'MlM,, (27) 
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with coefficients c""^. For 77 (9(1), c~ C'(l) by NDA. This will induce mixings between different 
components of the mesons. 

Consider one example K = (A^/M^)MqM'. By a field redefinition of Mq that diagonalizes 
the kinetic terms, one would have in the superpotential an 5*0(5) symmetry breaking operator 



where the coefficient c' ~ 0{1/Att) by NDA. Even though M' is not charged under the R symme- 
try, its scalar mixes with the i?-charge two singlet S' through the mass terms in the potential as 
A?M'S' . Thus the M' scalar can decay into two /2-axions while its fermionic partner can decay 
into an i?-axion and a gravitino. The estimate of the lifetime is similar to that of the DM in 6-1 
model. Another possibility for the scalar DM is that M' — )■ 5*0 5*0 if it is kinematically allowed, 
mass(M') > 2 mass(5'o). Here the Sq scalars in the final state will subsequently decay to i?-axions. 

5 Gauge mediation 

Having shown that viable DM can exist in the DSB sector, we now briefly discuss how to incor- 
porate such models into realistic gauge mediation scenarios. The goal is to introduce messenger 
fields with SM gauge charges which directly feel the SUSY breaking of the DSB sector. Yukawa 
couplings of messengers to the composite fields in the DSB sector are typically highly suppressed 
since such couplings are high dimension operators from the point of view of the microscopic DSB 
fields. Thus, we are led to introduce messengers which are directly charged under the hidden 
gauge group. 

For instance, in the 6-1 model from Sec. [3], one could add messengers charged under the strong 
SU{'o) group [111 113] III Specifically, we add messengers P, P which transform as a fundamental 
and anti- fundamental under SU{Q) with a superpotential mass term mpPP. For a SUSY breaking 
scale y/F ~ 10 — 100 TeV, mp must taken to be close to the strong-coupling scale A. At the 
scale A, one expects the SU{6) strong dynamics to generate Kahler operators of the form |Pp|$p 

^One could also try adding a pair of messengers charged under the gauged U{1) in the 6-1 model. However, 
since U (1) is weakly coupled at the DSB scale, the mass splittings of such messengers would be one-loop suppressed 
from A. Consequently the MSSM soft masses are two-loop suppressed from A, in essence yielding a two-step weakly 
coupled gauge mediation. To have acceptable soft masses, A has to be (10^ — 10^) TeV, which is unacceptably 
large from the point of the view of achieving DM in the DSB sector. 



6W = —S,M'. 



(28) 



This forces the fields charged under the global symmetry to obtain tiny VEVs 




(29) 
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where $ are the composite moduh. Such operators lead to SUSY breaking mass sphttings inside 
the messenger sector of order Am ^ 0{A). 

These messengers transmit SUSY breaking from the DSB sector to the SSM because they 
are charged under the SM gauge group, e.g. P + P transform as 5 + 5 under the GUT group 
SU{5)sM- With this messenger content, SU{6) is still asymptotically free, so the desired DSB 
dynamics persists. The low energy MSSM soft masses can be parametrized by the current-current 
correlators of the hidden sector as in general gauge mediation |13]. By NDA, 



where we used F ~ A'^/Att and C2{r) is the quadratic Casimir invariant of the representation r 
of the scalar. Here the factors of Att arise from NDA and are coincident with the loop factor in 
the weakly-coupled case. A ~ (10 - 100) TeV gives the MSSM gaugino mass around (100 GeV - 
1 TeV), as required. 

Note that adding these six sets of 5 + 5 messengers will ruin perturbative GUT unification. 
Of course, preservation of perturbative grand unification is not a necessary requirement, and in 
fact, DSB sectors with strongly coupled messengers tend not to have this property. However, 
in the context of the decaying DM scenarios envisioned in Sec. 12. 3[ perturbative unification is a 
desirable feature. The reason is that to get DM with cosmological lifetimes, the higher dimensional 
symmetry breaking operators must be suppressed by the GUT scale, and without perturbative 
unification, the motivation for such operators is lost. Possible ways of preserving perturbative 
gauge unification despite the large number of messenger fields have been discussed in Ref. |15] . 

One marginal possibility in the class of dynamical superpotential models is to consider the 
5-2 construction, i.e. the DSB model SU{n) x SU{2) x f/(l) (n odd) with = 5. Such a model 
can also yield strongly coupled DM, but need only have five 5 + 5 messengers to achieve gauge 
mediation. With A ~ (10 — 100) TeV, unification is marginally preserved |16]. We emphasize, 
though, that having DM decay is not necessary from the point of view of having a viable DM 
scenario, but is simply an interesting possibility in light of indirect DM detection. 

Turning to models with a quantum moduli space, note that adding 5 + 5 messengers does 
not work for the simplest ITIY model in Sec. HI as it ruins the asymptotic freedom of the gauge 
group SU{2). To add strongly coupled messengers, we must consider an ITIY- type model with 
a larger gauge group, for example an ITIY model based on a Sp{2N) gauge group with 2A^ + 2 
fundamentals Q and A^(2A^— 1) gauge singlets [H]. The effective superpotential is a generalization 
of the one for the SU{2) model (see Eq. (35) in Ref. [H]) and can be analyzed in the same way. For 

> 1, asymptotic freedom is preserved after adding messengers P{P) in the (anti)fundamental 




(30) 
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representation of Sp{2N) ^2]. = 2 is the unique choice for which perturbative unificiation is 
preserved after adding the messengers. 

Finally, messenger fields may themselves influence cosmology. The simplest messenger sectors 
preserve a messenger parity, so the lightest component is stable and may provide an additional 
DM component. Since these messengers are charged under the DSB gauge group, they experience 
strong annihilation cross sections, and as long as the lightest messenger is heavier than some 
unstable fields in the DSB sector, their relic abundance will be comparable to the other stable 
DSB fields. However, since the messengers are charged under the SM gauge group, there will 
always be a colored component, and there are strong constraints on colored particles with lifetimes 
r > 10^^ sec as they form exotic atoms [17]. Breaking the messenger parity with renormalizable 
interactions between the messenger and observable sectors could reintroduce the flavor problem. 
An option which preserves messenger parity is to make the triplet inside 5 + 5 decay through a 
dimension-five operator suppressed by the Planck scale, e.g., / (f9P^5^, where 5 refers to a SM 
multiplet t22j. 

6 Conclusions 

Most studies of DM focus on the weakly-interacting DM paradigm, where DM has mass around 
the electroweak scale ~ 0(100 GeV — 1 TeV) and couplings of weak interaction strength 

K ~ 0{1), leading to an annihilation cross section that yields the desired thermal relic density. 
This paradigm, however, is difficult to incorporate into another theoretically attractive paradigm, 
GMSB, which solves the flavor problem of the SSM. In standard GMSB, the LSP is the super-light 
gravitino, which cannot play the role of cold DM in the universe. 

In this paper, we have explored an alternative mechanism via which GMSB could still yield 
thermal relic dark matter, by considering DM in the DSB sector. Since the relic density of DM is 
controlled by the DM annihilation cross section, and since this cross section depends mainly on 
the ratio K^/mDM, it is possible to arrange uibm to be related to a heavier scale than the weak 
scale provided the DM couplings are correspondingly increased. Unitarity of the annihilation 
cross section then bounds the DM mass to be (9(10 — 100 TeV) and couplings to be of order a 
few to An. A DSB sector with the SUSY breaking scale a/F ~ (9(10 — 100 TeV) naturally realizes 
this situation, and heavy DM could be a composite state of SUSY breaking strong dynamics. 

We have found that DM can indeed arise from two of the most popular DSB scenarios: chiral 
models with a dynamically generated superpotential and vector-like models with a quantum 
moduli space. Both of these setups allow the DSB sector to have a large global symmetry that 
is preserved by SUSY breaking, and the lightest state(s) charged under this global symmetry 
are stable DM candidates. This setup also provides an additional probe of the dark sector: the 
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i?-axion. DSB models often preserve an approximate R symmetry which is broken when SUSY is 
broken. The resulting i?-axion is naturally light, 1 MeV - 10 GeV, with an irreducible contribution 
to its mass from supergravity effects. We have also argued that these DSB sectors can incorporate 
strongly-coupled messenger fields as necessary to have a complete model of GMSB. 

These models are an explicit realization of the scenario envisioned in Ref. [10], where decay 
of heavy composite DM might be responsible for various astrophysical anomalies. Indeed, the 
symmetries that stabilize DM need not be respected by high scale physics at the unification or the 
Planck scale. The compositeness of DM then forces the leading symmetry- violating operators to 
be high-dimensional, and consequently DM decays on cosmological time scales. In particular, we 
have found explicit DSB models where the leading symmetry violating terms are dimension-six 
operators. When suppressed by the Planck scale, such operators yield a DM lifetime of 10^^ sec 
for rriDu ~ 10 TeV, a range that can be probed by current and future DM indirect detection 
experiments. Since the leading DM decay is to /2-axions, and since by kinematics the i?-axion 
decays dominantly to SM leptons and pions without producing many hadrons, such a scenario 
can explain the lepton-rich astrophysical anomalies. 

The possibility of DM in the DSB sector presents a new direction for SUSY model building. 
While the bulk of previous DSB studies have focused on simply establishing SUSY breaking — the 
zeroth-order requirement for a DSB sector — it is worthwhile to develop tools to analyze the sym- 
metries and spectra of the SUSY breaking vacuum, especially for strongly-coupled models. We 
have shown that explicit DSB model have residual global symmetries and viable cold DM candi- 
dates, and we expect this to be true for a wide variety of DSB scenarios. In our models, besides 
the SSM spectrum which encodes limited DSB information, light states such as i?-axions provide 
additional handles on hidden sector dynamics which could be tested in near future astrophysics 
and collider experiments [TOl |26l B8] . 
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